home *** CD-ROM | disk | FTP | other *** search
/ Graphics Plus / Graphics Plus.iso / general / fractal / kaos.lha / eigenlib / cdiv.c < prev    next >
Encoding:
Text File  |  1989-11-18  |  304 b   |  19 lines
  1. cdiv(ar,ai,br,bi,cr,ci)
  2.  
  3. double ar,ai,br,bi,*cr,*ci;
  4. {
  5.  
  6. /*    complex division, (cr,ci) = (ar,ai)/(br,bi)   */
  7.  
  8. double s,ars,ais,brs,bis,fabs();
  9.  
  10. s = fabs(br) + fabs(bi);
  11. ars = ar/s;
  12. ais = ai/s;
  13. brs = br/s;
  14. bis = bi/s;
  15. s = brs*brs + bis*bis;
  16. *cr = (ars*brs + ais*bis)/s;
  17. *ci = (ais*brs - ars*bis)/s;
  18. }
  19.